Transport of Orbital- Angular-Momentum Entanglement through a Turbulent 

Atmosphere 
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We demonstrate experimentally how orbital-angular-momentum entanglement of two photons 
evolves under influence of atmospheric turbulence. We find that the quantum channel capacity is 
surprisingly robust: Its typical horizontal decay distance is of the order of 2 kilometers, demonstrat- 
ing the potential of photonic orbital angular momentum for free-space quantum communication in 
a metropolitan environment. 
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Quantum communication by means of entangled pho- 
ton pairs allows for an intrinsically secure transmission 
of data, by distributing the pairs via a free-space or fiber 
channel to distant parties Most popular is polariza- 
tion entanglement, which has dimensionality 2. Higher 
dimensionalities can be achieved using orbital-angular- 
momentum (OAM) entanglement @, or energy-time 
entanglement ; this route provides for a larger chan- 
nel capacity and an increased security against eavesdrop- 
pers However, the performance of a practical high- 
dimensional quantum channel is an open issue. Here, we 
address this issue for the case of OAM entanglement dis- 
tribution via a free-space channel. 

For quantum communication to be of practical rele- 
vance, it is imperative that the entanglement between 
the photons carrying the information survives over a rea- 
sonably long propagation distance. Entanglement distri- 
bution over fiber-based transmission lines has proven to 
maintain coherence over tens of kilometers [ij [l^, 11 1. 



However, the use of free-space links cannot be obviated 
when considering such purposes as airplane and satel- 
lite quantum links or hand-held communication devices 

3 51 0. 

The increased quantum-channel capacity that is avail- 
able when encoding the information in the OAM of 
the entangled photons was anticipated to be severely 
limited in a practical free-space link, due to atmo- 
spheric turbulence that causes wave front distortions. 
Several theoretical studies have addressed this aspect 
[3 EE E3, 0, H, 0, EH, but there is no unanimity 
on exactly how sensitive OAM entanglement is to atmo- 
spheric perturbations. So far, no experimental verdict 
has been obtained to clarify this issue. 

In this Letter, we present the first such experiment. We 
start with bipartite OAM entanglement of dimensional- 
ity 6, and demonstrate how the corresponding quantum 
correlations evolve when one of the photons traverses a 
turbulent atmosphere, emulated by controlled mixing of 
cold and hot air. Our experimental results are in ex- 
cellent agreement with our theoretical model, which is 



based on a Kolmogorov description of atmospheric tur- 
bulence. Specifically, we show how increasing strengths 
of turbulence degrade the Shannon dimensionality, which 
was introduced in Ref. [4] as a simple measure to quan- 
tify the channel capacity. Scaling up our system to real- 
life dimensions, we find that its typical horizontal decay 
distance is about 2 km. 

Our experimental setup is depicted in Fig. [1] The 
PPKTP crystal emits correlated photon pairs with 
complementary OA M, in a state of the form |\&) = 
J2i, p (i\hp) \~hp) [H, H3- Here, \l,p) indicates the 
Schmidt mode containing one photon with orbital an- 
gular momentum lh, with p the radial mode index, and 
we can write (r\l,p) = (r\l,p)e /y2w. The total num- 
ber of entangled azimuthal modes, the so-called angular 
Schmidt number K, is of order 30 The correlated 

photons are spatially separated by a 50/50 beam splitter. 

The entanglement is analyzed by means of two state 
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FIG. 1: Experimental setup. A type-I PPKTP crystal emits 
two frequency-degenerate photons (A = 826 nm) that are en- 
tangled in their OAM degree of freedom. A beam splitter 
serves to separate the twin photons spatially. The entangle- 
ment is analyzed by two angular-phase-plate projectors, vari- 
ably oriented at a and j3, respectively, which are linked to a 
coincidence circuit. Each phase plate has one elevated quad- 
rant sector with optical thickness A/2 (inset). In one of the 
beam lines we place a turbulence cell. 
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projectors, which are composed of an angular phase plate 
that is lens-coupled to a single-mode fiber and a single- 
photon counter. The two phase plates are identical and 
carry a purely azimuthal variation of their optical thick- 
ness: they have one elevated quadrant sector with an 
optical thickness that is A/2 larger than that of the re- 
mainder of the plate (see inset Fig. [T]) . The two phase 
plates can be rotated around their normals over an an- 
gle a and (3, respectively. The detection state |A(a)) (or 
|.B (/?))) of one such analyzer is a high-dimensional su- 
perposition of OAM modes, the relative phases of which 
depend on the orientation of the plate. It can be writ- 
ten as (r\A(a)) = {2/w )exp(~r 2 /w 2 )J2 l V^ie a(9+a \ 
where the Gaussian factor describes the fiber mode pro- 
file with field radius wq, and the summation over the 
orbital-angular-momentum states describes the phase im- 
print imparted by the phase plate. When rotating the 
quadrant phase plate over 27r rad, the analyzer scans a 
mode space of dimensionality D 
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In one of the beam lines, we place a turbulence cell 
where cold and hot air are mixed to bring about ran- 
dom variations of the refractive index that vary over time 
(Fig. [2a)). We can tune the strength of the turbulence 
by varying the heating power and air flow through the 
cell. Similar cells have been used as a realistic emulation 
of atmospheric turbulence [26| . Figure Wljo) gives an im- 
pression of the cell's functioning: We inject one of the 
analyzers backwards with diode laser light and monitor 
the beam, which traverses the turbulence cell, in the far 
field. We do this for two cases; the analyzer is equipped 
with no phase plate (top row), or with the quadrant- 
sector phase plate (bottom row). We observe that the 
input beams (left column) become deformed by the re- 
fractive index fluctuations, as can be seen when taking a 
10 ms snapshot (middle column). Time averaging these 
fluctuations over 10 s reveals a beam broadening that is 
spatially isotropic (right column). 

We describe our cell by the Kolmogorov theory of tur- 
bulence 27|. This standard model treats the optical 



effects of the atmosphere at any moment as a random 
phase operation e l ^ r \ the time evolution of which fol- 
lows a Gaussian distribution. It is conveniently described 
in terms of its coherence function, given by 



,i0(ri)— i0(r2) 



46.88 



(1) 



where (■ ■ -)t denotes averaging over time [28|. The rele- 
vant parameter in this model is the Fried parameter ro, 
being the transverse distance over which the beam pro- 
file gets distorted by approximately 1 rad of root-mean- 
square phase aberration [28|. In the absence of turbu- 
lence ro — *■ oo, but when turbulence becomes stronger, 
the spatial coherence is reduced and hence r shortens. 
From the Gaussian beam broadening in Fig. HJb) (top 
row) we can determine the relation between the Fried 
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FIG. 2: Beam corruption after passage through the turbu- 
lence cell, (a) Impression of how an OAM eigenmode, having 
a helical wavefront, gets distorted when transiting the tur- 
bulence cell. The cell consists of a 7 cm long, 26 mm diam- 
eter glass tube, containing several resistors that produce up 
to 60 W of heat. A gentle flow of room temperature air is 
driven through the tube, (b) Far-field intensity patterns of 
the analyzer, which is fed backwards with diode laser light at 
826 nm. The analyzer is equipped with no phase plate (top 
row), or quadrant-sector phase plate (with its sector aligned 
along the Cartesian axes) (bottom row). The diffraction lim- 
ited patterns (left column) get perturbed when turbulence is 
switched on (middle column): for mild turbulence, the dom- 
inant effect is a randomly evolving beam deflection; for the 
more severe turbulence conditions used here (wo/ro — 0.65), 
the beam profile can get significantly distorted. Taking a 10 
s time average reveals an isotropic beam broadening (right 
column). The apparent asymmetry along the diagonal in the 
bottom left and right windows is due to the 3% discrepancy 
of the quadrant phase step from the ideal value of tt. 



parameter and the 1/e beam size Wq, 

WQ _ \f{wie/w d l) 2 - 

r ~ 3.0 



(2) 



with Wdi and wi e the 1/e far-field radius of the diffrac- 
tion limited beam and long-exposure broadened beam, 
respectively |29|. 

We calculated the effect of Kolmogorov turbulence on 
our analyzer state |A(a)). The blue curve in Fig. [3]^a) 
shows the survival probability of an OAM eigenmode 
I = Iq upon passing through a turbulent atmosphere 
as described by Eq. ([1]). The survival probability de- 
grades gradually for increasing turbulence strength. We 
note that this decay depends on the ratio u>o/ro only 
and not on the specific OAM eigenvalue Iq, provided 
that the propagation distance L is small compared to the 
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FIG. 3: Mode scattering due to turbulence, (a) Survival prob- 
ability of an analyzer's OAM eigenmode / = lo as a function 
of turbulence strength (blue). The red and green curves de- 
note turbulence-induced coupling probabilities to neighboring 
modes for Al = ±1 and Al — ±2, respectively, (b) Spreading 
of the lo OAM eigenmode (blue bar) over its neighbors for 
wo/ro = 0.65 (red bars). 



diffraction length zr — ttwq/X. Furthermore, the turbu- 
lence produces a coupling between the orthogonal OAM 
modes, leading to a non- vanishing mode overlap between 
the Iq eigenmode and its neighbors Al = ±1 (red) and 
Al = ±2 (green). A different perspective on this mode 
mixing is presented in histogram Fig. [3jb), which shows 
how an OAM eigenmode (blue bar) spreads out over 
its neighboring azimuthal modes for wo/ro = 0.65 (red 
bars). We note that normalization is not preserved, be- 
cause some intensity is scattered to radial modes that 
are not sustained by the single-mode fiber. This illus- 
trates the importance of taking into account the radial 
content of the generated two-photon state and the ana- 
lyzers' detection states when dealing with OAM modes 
in the presence of turbulence. 

In the experiment, the phase plates are rotated around 
their normals, and the photon coincidence probability 



P(a-/3) = (\(A(a)\(B((3)\S A m' 



(3) 



is recorded as a function of their independent orienta- 
tions. Here, the time-averaged behaviour of the turbu- 
lence scattering operator Sa is known by means of its 

coherence function Eq. (T]): /s A \A(a)){A(a)\S A \ = 

/dr 1 dr 2 |r 1 )(r 1 |A(a))(A(a)|r 2 )(r 2 | ( e «(rO-*(*»)) t . Fig _ 
ure 0] shows our main experimental results. In the ab- 
sence of turbulence, we observe a piecewise-parabolic co- 
incidence curve (blue circles), i.e. the coincidence rate 
follows a parabolic dependence for \a — /3\ < ir/2 and 
is zero elsewhere The coincidence rate depends on 
the relative orientation of the phase plates only. We 
have investigated how the coincidence rates evolve for 
6 turbulence strengths, two of them shown in Fig. SJ 
wo/ro = 0.30 (green triangles) and Wo/ro = 0.65 (red 
stars). The latter strength was also used for Fig. [2jb) 
and^b). Note that the 20 s integration time assures 
isotropic sampling of the wavefront fluctuations (see Fig. 
H^b)). We observe a partial "smoothening" of the coin- 
cidence curve, which is excellently described by our the- 




FIG. 4: Survival of OAM entanglement under influence of 
turbulence. Experimental coincidence rates (data points) and 
theoretical predictions (curves) obtained with two quadrant- 
sector phase plates for: no turbulence (blue), wo/ro = 0.30 
(green) and wo/ro = 0.65 (red). The inset shows a blow-up 
of the wiggles around a — j3 — it/2. 



oretical predictions based on Eqs. ([J) and (J3j) , without 
any fit parameter. The turbulence-induced wiggles at 
| a — j3\ = 7r/2 are reproduced remarkably well (see in- 
set). 

Figure 2] shows that the coherence of the entangled 
state is partly conserved. We quantify its robustness 
in terms of the Shannon dimensionality D, which is 
a measure of the quantum channel capacity [4j. It is 
an operationally defined measure and gives the effec- 
tive number of modes the combined analyzers have ac- 
cess to when scanning over their possible settings, viz. 
the phase-plate orientations. For two identical analyz- 
ers in the absence of turbulence, we can express D in 
terms of their pure detection state operators pa and ps, 
where p A = \A(a))(A(a)\, as D = l/Tr({p A ) a (pb)p) = 
l/Tr((pA) a ) 2 - Here, (pa) a is the density matrix averaged 
over all phase plate orientations a. 

In the presence of turbulent scattering, however, the 
detection state becomes randomly time dependent: pa — 
PA(t)- The relevant detection state operator is therefore 
not pa, but rather {pa)i, i-e., the density matrix aver- 
aged over time. In general, {pa)i is no longer a single- 
mode projector, but a mixed positive operator. In other 
words, when averaging over the random fluctuations, the 
detection state becomes multimode, which degrades the 
analyzer's modal resolution. Therefore, in the presence 
of turbulence, the Shannon dimensionality for mixed de- 
tection states is given by 



D = 



Tr((p)t) 2 
Tr((p)t,a) 2 



(4) 



where (p)t,a denotes the average of the detection state 
operator p over time t and orientation a. The inverse of 
the numerator of Eq. (0| can be interpreted as the effec- 
tive number of modes captured by the analyzer for fixed 
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FIG. 5: Decay of Shannon dimensionality. Experimental 
(data points) and theoretical (curves) dimensionality as a 
function of turbulence strength, for two quadrant-sector phase 
plates (circles) and two half-sector phase plates (triangles). 
The turbulence strength is expressed in the ratio wo/ro (lower 
abscissa) and in the corresponding real-life propagation dis- 
tance L (upper abscissa). 



orientation a. Alternatively, it also represents the purity 
P = Tr((p) t ) 2 < 1 of the detection state, whose evo- 
lution follows the survival probability discussed in Fig. 
[3Ja) and is independent of the specific phase plates in 
use. Moreover, the numerator gives the maximum coin- 
cidence rates obtained at the peaks in Fig. [4] up to a 
constant scaling factor [30] . 

Experimentally, D can straightforwardly be extracted 
from the coincidence curves in Fig. 2) D = 2wN max /A, 
where N max is the maximum coincidence rate and A is 
the area underneath the curve. Figure [5] shows how D 
evolves for increasing turbulence strength according to 
theory and experiment. In the absence of turbulence, 
we find an experimental value D — 5.7 vs. a theoretical 
prediction D = 6 4]. As the turbulence strength in- 
creases, the modal resolution of the analyzers degrades, 
constraining the dimensionality to smaller values, ulti- 
mately to D = 1. The number of communication modes 
is reduced by ~ 50% to D = 3.1 when wq/tq = 0.65 
(theory D — 3.3). Considering the severity of the wave 
front distortions (see Fig. HJb)), we conclude that the 
channel capacity is surprisingly robust. For comparison, 
we also plotted our experimental results obtained with 
two half-sector phase plates, having one semicircle phase 
shifted by it (see inset in Fig. [5]). For this case we ob- 
serve that the dimensionality, initially at a value D = 3, 
decays considerably more slowly. This indicates that the 
resilience to atmospheric turbulence is very sensitive to 
the nature of the OAM superposition state, an aspect 
also noted in Ref. 19]. We expect that the search for 
such optimal states will lead to even more tenacity of the 
channel capacity. 



The results presented here allow us to estimate the 
horizontal propagation distance L that can be reached 
for free-space quantum communication outside the lab- 
oratory, since the Kolmogorov theory (Eq. (T])) used to 
describe our data is also a fair description of a real-life 
atmosphere [27J. Let us consider the 50% decay level 
of the channel capacity, roughly reached in our exper- 
iment for wo/ro — 0.65 (D — 3.1). For horizontal 
propagation, the Fried parameter can be expressed as 
ro = 3.02(fc 2 LC 2 ) -3 / 5 [3lJ. Assuming moderate ground- 
level perturbations (C 2 = 10 -14 m~ 2 / 3 ) [32j, a wave- 
length A = 1550 nm in the transmission window of the 
atmosphere and a beam size wo — 6 cm, we find a 
propagation length of 2 km (satisfying the requirement 
L < Zr). This is sufficient for quantum key distribu- 
tion in a metropolitan environment, offering the advan- 
tages of high-dimensional entanglement as compared to 
the case of 2D polarization entanglement. We note that 
this distance could be enhanced significantly if one were 
to incorporate additional adaptive optics techniques [33| . 
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